GLOBAL WEAK SOLUTION AND LARGE-TIME BEHAVIOR FOR 
THE COMPRESSIBLE FLOW OF LIQUID CRYSTALS 



DEHUA WANG AND CHENG YU 

Abstract. The three-dimensional equations for the compressible flow of liquid crystals 
are considered. An initial-boundary value problem is studied in a bounded domain 
with large data. The existence and large-time behavior of a global weak solution are 
established through a three-level approximation, energy estimates, and weak convergence 
for the adiabatic exponent 7 > |. 



1. Introduction 

In this paper, we consider the fohowing hydrodynamic system of partial differential 
equations for the three-dimensional flow of nematic liquid crystals ([^ [TH [T7]): 

pt + d\v{pu) = 0, (1.1a) 
{pu)t + div(pu u) + VP{p) = pAu - Adiv (vd Vd - (^|Vdp + F(d))/3^ , (1.1b) 
dt + u-Vd = e(Ad-/(d)), (1.1c) 

where p > denotes the density, u G the velocity, d G the direction field for the 
averaged macroscopic molecular orientations, and P = ap^ is the pressure with constants 
a > and 7 > 1. The positive constants p,X,6 denote the viscosity, the competition 
between kinetic energy and potential energy, and the microscopic elastic relation time for 
the molecular orientation field, respectively. The symbol denotes the Kronecker tensor 
product, I3 is the 3x3 identity matrix, and Vd Vd denotes the 3x3 matrix whose 
ij-th entry is < 5^.d,c?a;^d >. Indeed, 

Vd0 Vd = (Vd)TVd, 

where (Vd)""" denotes the transpose of the 3x3 matrix Vd. The vector-valued smooth 
function /(d) denotes the penalty function and has the following form: 

/(d) = VdF(d), 

where the scalar function -F(d) is the bulk part of the elastic energy. A typical example is 
to choose F{d) as the the Ginzburg-Landau penalization thus yielding the penalty function 
/(d) as: 

F{d) = -L(|d|2 - 1)2, /(d) = Jl^ddp - l)d, 
4(7(3 2cJo 
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where ctq > is a constant. We refer the readers to [21 HI [5l [TTl [161 E] ^or more physical 
background and discussion of hquid crystals and mathematical models. 

There have been many mathematical studies on the incompressible flows of liquid crys- 
tals. In Lin-Liu |17l \TE[ \T9\ I20j . the global existence of weak solutions with large initial 
data was proved under the condition that the orientational configuration d(x, t) belongs to 
H^, and the global existence of classical solutions was also obtained if the coefficient fi is 
large enough in the three-dimensional spaces. The similar results were obtained also in [26] 
for a different but similar model. The global strong solution was established in Hu-Wang 
|14j . When the weak solutions are discussed, the partial regularity of the weak solution 
similar to the classical theorem of Caffarelli-Kohn-Nirenberg [1] was obtained in |20J (and 
also [12]). The existence of weak solutions to the density-dependent incompressible flow 
of liquid crystals was proved in fl5j . The compressible flow (jl.ip of liquid crystals is much 
more complicated and difficult to study mathematically due to the compressibility. In the 
one-dimensional case the global existence of smooth and weak solutions to the compressible 
flow of liquid crystals was obtained in [3]. Our aim of this paper is to establish the global 
existence of weak solutions {p, u, d) to the three-dimensional compressible flow (jl.ip of 
liquid crystals in a bounded smooth domain C i?^, with the following initial-boundary 
conditions: 

{p,pu,d)\t=o = {po{x),mo{x),do{x)), x e n, (1.2) 

and 

u\dn = 0, d\gn = do(x), x G dQ, (1.3) 

where 

PoeL''{n), po>0; do€ L^{n)nH\n); 

mo€L^(n), mo = Oif/9o = 0; e L^(J7). 

Po 

When the direction field d does not appear, (jl.ip reduces to the compressible Navier- 
Stokes equations. For the compressible Navier-Stokes equations. Lions in [22] introduced 
the concept of renormalized solutions to overcome the difficulties of large oscillations and 
proved the global existence of finite energy weak solutions for 7 > 9/5, and then Feireisl, et 
al, in [SlITllH] extended the existence results to 7 > 3/2. Hu-Wang in [13] adopted Feireisl's 
techniques to obtain global existence and large-time behavior of weak solutions with large 
initial data for the magnetohydrodynamics. In this paper we shall study the initial- 
boundary value problem ()l.ip - (ll.3p for liquid crystals and establish the global existence 
and large-time behavior of weak solutions for large initial data in certain functional spaces 
with 7 > 3/2. To achieve our goal, we will use a three-level approximation scheme similar 
to that in [6l[7j, which consists of Faedo-Galerkin approximation, artificial viscosity, and 
artificial pressure. Then, motivated by the work of [7], we will show that the uniform 
bound of the density p'^'^" in for some a > ensures the vanishing of artificial pressure 
and the strong compactness of the density. To overcome the difficulty of possible large 
oscillation of the density, we adopt the idea of Lions and Feireisl in [H El [19] based on 
the weak continuity of the effective viscous fiux for the Navier-Stokes equations. For our 
equations (jl.ip of liquid crystals, the effective viscous flux is P — ^divu. For this purpose, 
we also need to develop some estimates to deal with the direction field and its coupling 
and interaction with the fluid variables. To deal with the equation (jl.lcp for a given u, we 
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will follow the same idea of Hu-Wang [13j to establish the solvability of the direction field. 
It is crucial to obtain sufficiently strong estimates on the direction field d to recover the 
original system (II. Ih . We will derive an energy inequality from (jl.ip directly, but it can 
not provide us with sufficient regularity for the direction field d. Thanks to the Gagliardo- 
Nirenberg inequality and the maximum principle, we deduce that Vd G -^^^((0, T) x 17) for 
any T > which can be used to control the strongly nonlinear terms containing Vd in 
Then we finally establish the existence of global weak solution to (|l.ip - (jl.3p . We 
just noticed that a similar existence result was obtained independently in [25j. Motivated 
by [8j and [13] , we will also establish the large-time behavior of the global weak solutions. 

We organize the rest of the paper as follows. In Section 2, we deduce a prior estimates 
from (jl.ip . give the definition of finite energy weak solutions, and also state our main 
results. In Section 3, we discuss the solvability of the direction vector d in terms of 
u. In Section 4, we establish the global existence of solutions to the Faedo-Galerkin 
approximation to (jl.ip . In Section 5 and Section 6, we use the uniform estimates to recover 
the original system by vanishing the artificial viscosity and artificial pressure respectively. 
In Section 7, we prove the large-time behavior of the global weak solutions. 



In this section, we derive some basic energy estimates for the initial-boundary problem 
(|l.ip - (11.3p . introduce the notion of finite energy weak solutions in the spirit of Feireisl 
[Sl E], and state the main results. 

Without loss of generality, we take 9 = a = 1. First we formally derive the energy 
equality and some a priori estimates, which will play a very important role in our paper. 
Multiplying ()l.lbp by u, integrating over Q, and using the boundary condition (jl.3p . we 
obtain 



2. Energy Estimates and Main Results 




Using the equality 



div(Vd0 Vd) = V(-|Vd 



) + (Vd)' - Ad, 



we have 




Hence, we obtain 




(2.1) 
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Multiplying by A(Ad — /(d)) the both sides of (jl.lcp and integrating over Q, we get 
-dt [ X^-^^dx -dt [ XF{d)dx - [ XVdF{d)udx + X [ (Vd)"^ • Ad • udx 

= X [ |Ad- /(d)|2dx. 
Jn 

Then, from ()2.ip . we have the following energy equality to the system (jl.ip . 

dt ^ Q/.|u|2 + + ^|Vd|2 + AF(d)) dx 

+ [ (/x|Vu|2dx + A|Ad- /(d)|2)dx ^^-^^ 
Jn 

= 0. 

Set 

and assume that E{0) < oo. From (|2.2p . we have the following a priori estimates: 

p\u\^GL°-i[0,T];L\m 

peL^{[0,T]-L''{n)); 
VdG L~([0,r];L2(J7)); 
F(d)GL-([0,r];Li(J^)); 
VuGL2Qo,r];L2(J7)); 

and also 

Ad-f{d)GL\[0,T];L\n)). (2.3) 

Although the above estimates will play very important roles in proving of our main exis- 
tence theorem, they cannot provide sufficient regularity for the direction field d to control 
the strongly nonlinear terms containing Vd. To overcome this difficulty, we need the 
following lemma (see [9]): 

Lemma 2.1. If there exists a constant Cq > such that d-/(d) > for all \d\ > Cq > 0, 
then d G L°°((0,T) x Q), Vd G L^{{0,T) x Q). 

Proof. On one hand, if |d| < Co, we have 

d G L°°((o,r) X n). 

On the other hand, if |d| > Cq, taking the scalar product of equation (jl.lcp with d yields 
^ildp - A|dp + u • V|d|2 + 2d • /(d) + 2|Vdp = 0, 

which implies 

at|d|2- A|dp + u- V|d|2 < 0. 
Using the maximum principle for |dp to obtain 

d G L°°((0,r) X f)). (2.4) 
Using ()2.3p . ()2.4p . smoothness of /, and together with elliptic estimate, we get 

d G L'^{0,T;H^{n)). 
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Using the Gagliardo-Nirenberg inequality, for some constant C > 0, 

1 i_ 

||Vd||i4(j^) < C||Ad|||2(f^)||d|||^(j^) + C7||d||i^(o), 

which means 

vd G L^((o,r) X n). 

The proof is complete. □ 

Through our paper, we will use C to denote a generic positive constant, V to denote 
Cq°, and P' to denote the sense of distributions. To introduce the finite energy weak 
solution (yO, u, d), we also need to take a differentiable function 6, and multiply (jl.lap by 
b' (p) to get the renormalized form: 

b{p)t + dw{b{p)u) + {b'{p)p - 6(p))divu = 0. (2.5) 

We define the finite energy weak solution (p, u, d) to the initial-boundary value problem 
(jl.ip - (|1.3p in the following sense: for any T > 0, 

•P>0, pe L^{[0,T];L^n)), u€L\[0,T];W^'\n)), 

de L°°((o,r) X n)nL°°{[o,T]-H\n))nL^{[o,T];H\n)), 

with pu, d)(0, x) = {po{x),mQ{x),dQ{x)) for x e fl; 

• The equations hold in D' {{0,T) x Q), and (fTTal) holds in D' {{0,T) x M^) 
provided p, u are prolonged to be zero on \ il; 

• ([23]) holds in D'{{0,T) x Q), for any b G C^{R+) such that 

6 (z) = for all z G large enough, say z > M, (2.6) 

where the constant M may vary for different function 6; 

• The energy inequality 

E{t)+ I I (p|Vup(ix + A|Ad-/(d)|2)dxds < -E(O) 

holds for almost every t G [0, T]. 

Remark 2.1. It's possible to deduce that (|2.5p will hold for any b G C^{0,oo) n C[0,oo) 
satisfying the following conditions 

\b'{z)\ < c{z°' + Z2 ) for all z > and a certain a G (0, ^) (2.7) 

provided (p, u, d) is a finite energy weak solution in the sense of the above definition (see 
details in [?])• 

Now, our main result on the existence of finite energy weak solutions reads as follows: 

Theorem 2.1. Assume that Q. is a hounded domain of the class C'^'^^ , z/ > 0, and 

7 > |. // there exists o Co > 0, such that d • /(d) > for all |d| > Co > 0. Then for any 
given T > 0, the initial-boundary conditions ()l.ip - ()1.3p has a finite weak energy solution 
(p,u,d) on (o,r) X n. 
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Remark 2.2. The typical example 

F(d) = -L(|d|2 - 1)2, /(d) = -L(|dp - i)d, 

4ao 2cJo 

satisfies the assumption of Theorem 12.11 with Cq = 1. Thus the theorem holds for the 
typical case. 

Motivated by [8j and [13], we establish the following result on the large-time behavior 
of the weak solutions to the problem (jl.ip - (|1.3p : 

Theorem 2.2. Assume that (/), u, d) is the finite energy weak solution to (|l.ip - (|1.3p given 
in Theorem \2.1l then there exists a stationary state of the density ps which is a positive 
constant, a stationary state of velocity Ug = 0, and a stationary state of direction field dg 
such that, as t —)• oo, 

p(t,x) — )■ Ps strongly in L'^($7), (2-8) 
u(t,x) Us = strongly in L'^{Q), (2-9) 
d{t,x) ds strongly in H^{n), (2.10) 



where dg solves the equation 
with the boundary condition 
and Ps satisfies the following relation: 



Ads = fids) (2.11) 
d.lef^ = do; (2.12) 



Vp] = -Adiv ( Vd, Vd, - (^|Vd,|2 + F(d,))/3 ) . (2.13) 



Remark 2.3. The existence and uniqueness of (j2.1ip and (j2.12p can be guaranteed from 
the elliptic theory, and d^ E C'^{Q) n C^{Q) by the standard elliptic estimates (see [TO]). 

Remark 2.4. Denote H = p2 — XF{ds), then equation (I2.13P can be rewritten as 

VH = -XVds ■ Ads. (2.14) 

By (j2.1ip and /(d^) = Vds-P'(ds), we can deduce that 

V{H + XF{ds)) = 0, 

that is, 

Vp] = 0. 

By the way, in general, we cannot solve the equation (I2.14p without the condition /(d^) = 
Vd.F(d.). 

Remark 2.5. Our asymptotic equations (j2.1ip - (l2.13p for the compressible flow of liquid 
crystals as t — )• oo are similar to those for the incompressible flow of liquid crystals obtained 
by Lin-Liu (see [i9\). In particular, the asymptotic equations (|2.1ip . (|2.12p . and (I2.14[) 
share the same form with those in |19j . 



GLOBAL SOLUTION TO THE FLOW OF LIQUID CRYSTALS 



7 



The proof of Theorem 12.11 is based on the following approximation scheme: 



pt + div(pu) = eAp, 

(pu)t + div(pu (g) u) + VP(p) + 5Vp^ + eVu • Vp 



(2.15a) 



= //Au - Adiv ^Vd Vd - (^|Vdp + F(d))/3 
dt + u - Vd = Ad - /(d), 



(2.15b) 



(2.15c) 



with appropriate initial-boundary conditions. Following the approach of Feireisl [6, 7J, we 
shall obtain the solution of (jl.ip when e — )• and 5 — )• in (j2.15p . We can solve equation 
(|2.15ap provided u is given. Indeed, we can obtain the existence by using classical theory 
of parabolic equation and overcome the difficulty of vacuum. Next we can also solve 
equation (j2.15cp when u is fixed. By a direct application of the Schauder fixed point 
theorem, we can establish the local existence of u, and then extend this local solution to 
the whole time interval. Note that the addition of the extra term eVu • Vp is necessary 
for keeping the energy conservation. The last step is to let e — )• and 5 — )• to recover 
the original system. We remark that the strongly nonlinear terms containing Vd can 
be controlled by the sufficiently strong estimate about Vd obtained from the Gagliardo- 
Nirenberg inequality. In order to control the possible oscillations of the density p, we 
adopt the methods in Lions |22] and Feireisl [6l [7] which is based on the celebrated weak 
continuity of the effective viscous fiux P — pdivu. We refer the readers to Lions [22], 
Feireisl [6n7j, and Hu-Wang ^13j for discussions on the effective viscous flux. 



To solve the approximation system (j2.15p by the Faedo-Galerkin method, we need to 
show that the following system can be uniquely solved in terms of u: 



which can be achieved by the two lemmas below. 

Lemma 3.1. If u G C([0, T]; Cq (Jl, M^)), then there exists at most one function 



which solves (13. ip in the weak sense on Q x (0,T), and satisfies the initial and boundary 
conditions in the sense of traces. 

Proof. Let di,d2 be two solutions of (|3.ip with the same data, then we have 



3. The Solvability of the Direction Vector 



dt + u • Vd = Ad - /(d) 
d|t=o = do, = do, 



(3.1a) 
(3.1b) 



d G L\0, T; Hl{n)) n L~([0, T]; H^ft)) 



(di - d2)t + u • V(di - da) = A(di - ds) - (/(di) - /(da)). 



(3.2) 
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Multiplying (j3.2p by A(di — d2), integrating it over 17, and using integration by parts and 
the Cauchy-Schwarz inequality, we obtain 

dt [ |V(di -d2)|2dx + 2 [ |A(di -d2)|2dx 
Jn Jn 

= 2 / (V(di - d2))^ • (A(di - d2)) ■udx + 2 [ (/(di) - /(d2))(A(di - d2))dx (3.3) 
Jn Jn 

<C I |V(di-d2)|2dx+ / |A(di-d2)|2dx, 
Jn Jn 

where we used the fact that / is smooth. Then 

dt I |V(di -d2)p(ix+ / |A(di-d2)|2dx<C / |V(di -d2)p(ix, (3.4) 
Jn Jn Jn 

and Lemma 13.11 follows from Gronwall's inequality, the above inequality, together with 
Lemma 12. li □ 

Lemma 3.2. Let C 6e a bounded domain of class C^"*"^, > 0. Assume that 
u G C([0,T]; Cq(0,M^)) is a given velocity field. Then the solution operator 

u ^ d[u] 

assigns to u £ C([0, T]; Cg (0; M'^)) the unique solution d of (j3.ip . Moreover, the operator 
u I — > d[u] maps bounded sets in C{[0,T]; Cq{CI;M.^)) into bounded subsets of 

Y := L\[0,T];H^in))nL°-i[0,T];H\n)), 

and the mapping 

is continuous on any bounded subsets of C{[0,T];Cq{Q;] 



Proof. The uniqueness of the solution to (|3.ip is a consequence of Lemma \'6.2\ and the 
existence of a solution can be guaranteed by the standard parabolic equation theory. 
By ()3.4p . we can conclude that the solution operator u i — > d(u) maps bounded sets in 
C{[0,T];Cq{Q;W^)) into bounded subsets of the set Y. Our next step is to show that 
the solution operator is continuous from any bounded subset of C([0, T]; Cq (Jl)) to Y. 
Let {un}'^^i be a bounded sequence in C{[0,T];C^{n)), that is to say, u„ G B{0,R) C 
C{[0,T];Cf{n)) for some R>0, and 

Un uin C{[0,T];C^{n)) as n ^ oo. 
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Here, we denote d[u] = d, and d[u„] = d„, so we have 
dt [ J|V(d„-d)|2dx+ / \A{dn-d)\^dx 

= / (u • Vd - u„ • Vd„)(V(d„ - d))dx + [ (/(d) - /(d„)) • (A(d„ - d))dx 

< /" (|u - u„| • |Vd| + |u„||V(d - d„)|)|A(d„ - d)\dx + C I |V(d„ - d)\^dx (3.5) 
Jn Jn 

<\\Un-u\\L^\\Vd\\l,+C\\V{d-dn)\\l2+l: [ \A{dn-d)\^dx 

^ Jn 

<C||u„-u||Loo + i||V(d-d„)||i2, 
where we used facts that d„ is bounded in Y and / is smooth. This imphes that 

\dt [ |V(d„ - d)\^dx + |A(d„ - d)fdx 

^ Jn Jn (3.6) 

< C\\Un-u\\Loo+C\\V{dn-d)\\l2. 

Integrating p.6p over time t £ (0, T), and then taking the upper hmit over n on the both 
sides, we get, noting that u„ — u in C([0,T]; Cq{^);. 



-Hmsup / |V(d„, - d)|^dx + - Hmsup /" [ \A{dn - d)\'^dxdt 
^ Jn ^ Jo Jn 

<Chmsup/ \\V{d„,-d)\\l2dt (3.7) 
" Jo 

fT 

<C hmsup ||V(d„ - d)||?2(it, 
Jo " 

thus, using Gronwall's inequahty to (j3.7p and noting that d„, d share the same initial 
data, we have 

hmsup / |V(d„ - d)pdx = 0, 
" Jn 

which means, from (|3.7p again, 

hmsup/ f \A{dn- d)\'^dxdt = 0. 
" Jo Jn 

Thus, we obtain 

d„ — > d in Y. 

This completes the proof of the continuity of the solution operator. □ 
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4. The Faedo-Galerkin Approximation Scheme 

In this section, we establish the existence of solution to the following approximation 
scheme: 

pt + div(pu) = eA/9, (4.1a) 
{p\i)t + div(pu (g) u) + VP(p) + 5Vp^ + eVu • Vp 

= /xAu - Adiv [vd Vd - (^|Vd|2 + F{d))I^ , (4.1b) 

dt + u-Vd = Ad-/(d), (4.1c) 

with boundary conditions 

V/9 • v\an = 0, (4.2a) 

dbn = do, (4.2b) 

ubn = 0, (4.2c) 

together with modified initial data 

p\t=Q = Po,5ix), (4.3a) 

pu\t=o = mo,5(x), (4.3b) 

d|t=o = do(x). (4.3c) 

Here the initial data pq^s{x) £ C'^(r2) satisfies the following conditions: 

0<S<po,six)<6~^^, (4.4) 

and 

Po,6{x) ^ pin L^{n), \{po,5 < poll ^0 as 5 ^ 0. (4.5) 

Moreover, 

mo,(x) = (™V/^7f)^ff)' (4.6) 
[0 if po,5{x) < Po{x). 

The density p = p[u\ is determined uniquely as the solution of the following Neumann 
initial-boundary value problem (see Lemmas 2.1 and 2.2 of [7]): 

Pt + div(/9u) = eAp, (4.7a) 

V/9-z.|aQ = 0, (4.7b) 

p\t=o = Po,5ix), (4.7c) 

To solve (I4.1bp by a modified Faedo-Galerkin method, we need to introduce the finite- 
dimensional space endowed with the Hilbert space structure: 

Xn = span{r]i)i^i, n E {1, 2, 3, • • • }, 
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where the Unearly independent functions rji £ 'D{Q)'^, i = 1,2, . . . , form a dense subset in 
Cq(J7,M^). The approximate solution u,„ should be given by the following form: 



(4.9) 



pun(r) • r/dx — / mo,5 • rjdx 
Jn 

- / / Adiv(Vd0Vd-(^|Vd|2 + F(d))/3)-r?dxdt 
Jo Jn 2 

for any t G [0,T] and any r] G X„, where £,5,(3 are fixed. Due to Lemmas 2.1 and 2.2 of 
[7] and our Lemmas 13.11 and 13.21 the problem (|3.ip . (j4.7|) and (j4.8|) can be solved at least 
on a short time interval (0, T„) with < T by a standard fixed point theorem on the 
Banach space C([0,T],X„). We refer the readers to [7J for more details. Thus we obtain 
a local solution (^„,u„,d„) in time. 

To obtain uniform bounds on u„, we derive an energy inequality similar to (12. 2p as 
follows. Taking r] = u„(t, x) with fixed t in (14. ip and repeating the procedure for a priori 
estimates in Section 2, we deduce a "Kinetic energy equality": 

dt f l\pnWn\^ + -^pl + -^pi + ^\Vdn\'^ + \F{dn)]dx + ^l f \Vvin\^dx 

+ \ f \/\dn- f{dn)\^dx + e I (7pr' + '^/3p^')|V/9„pdx = 0. 

Jn 

The uniform estimates obtained from (|4.9p furnish the possibility of repeating the above 
fixed point argument to extend the local solution u„ to the whole time interval [0,T]. 
Then, by the solvability of equation (14. 7p and ()3.ip . we obtain the functions (/?n,d„) on 
the whole time interval [0,r]. 

The next step in the proof of Theorem 12. 11 is to pass the limit as n — )• oo in the sequence 
of approximate solutions {yOn,u„,d„} obtained above. We observe that the terms related 
to u„ and p„ can be treated similarly to [7]. It remains to show the convergence of the 
terms related to d„. 

By (14. 9p . smoothness of /, and elliptic estimates, we conclude 

Vun(GL'{%T];L\^)), (4.10) 
Ad„ - /(d„) is bounded in ^^([0, T]-L^{a)), (4.11) 

and 

dn(^L°^{[Q,T]-Hlmr\L\[Q,T]-Hlm. 

This yields that 

Ad„ - /(d„) ^ Ad - /(d) weakly in L^{[Q,T]-L^{a)), 

and 

d„ ^ d weakly in L°°{[{),T];HI{9)) ^ L^{[Q,T];Hl{n)). (4.12) 
Using corollary 2.1 in [6] and (I4.1cp . we can improve (I4.12P as follows: 

d„^d inC([0,r];LLafc(f^))- 
Next we need to rely on the following Aubin-Lions compactness lemma (see [23]): 
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Lemma 4.1. Let Xq,X and Xi be three Banach spaces with Xq C X C Xi. Suppose that 
Xq is compactly embedded in X and that X is continuously embedded in Xi; Suppose also 
that Xq and Xi are reflexive spaces. For 1 < p,q < oo, let 

nil 

W = {ue LP{[0,T];Xo)\- G L^[0,T]; X,)}. 

Then the embedding ofW into LP{[0,T]; X) is also compact. 

We are now applying the Aubin-Lions lemma to obtain the 
From Lemma [2. H we have 

d„eL°°((0,T)xJ]), 

and 

Vd„ G L^((0,r) X 17). (4.13) 

Using (j4.1cp . we have 

Wdtdnh^n) < C||u„ • Vd„||i2(Q) + C||Ad„ - /(d„)||i2(f^) 

< C'||u„||^4(Q) + C||Vd„||^4(Q) 

+ C'l|Vd„|||4(Q) + C||Ad„ - /(d„)||i2(Q), 

where we used embedding inequality, the values of C are variant. Thus, (j4.10p . ()4.1ip and 
([TOD yield 

\\dt<in\\L'2{[0,T];L2{n)) < C. 

Notice that C C and the injection ^ is compact, applying Lemma |4. II 
we deduce that the sequence {d„}^]^ is precompact in L'^{[0,T]; (Q)). 

Summing up the previous results, by taking a subsequence if necessary, we can assume 
that: 

d„^d inC7([0,r];LLafc(O)), 

d„ ^ d weakly in L'^{0,T;H'^{n)) n L°^{0,T;H^{n)), 

dn^d strongly in L^{0,T; (n)), 

Vdn Vd weakly in L^((0,T) x 

Ad„ - /(d„) ^ Ad - /(d) weakly in L\o,T; L\n)), 

F{dn) F{d) strongly in L'^{0,T; H^{n)). 

Now, we consider the convergence of the terms related to d„ and Vd„. Let if he a test 
function, then 

/ (Vd„ Vd„ - Vd Vd) • V^pdxdt 
Jn 

< / (Vd„ Vd.„ - Vd„ Vd)V((Cdxdt + / (Vd„ Vd - Vd Vd)Vv?(ij;dt ^'^'^'^^ 
Jn Jn 

< C|| Vd„||i2(Q) ||Vd„ - Vd||i2(Q) + C|| Vd||i2(Q) ||Vd„ - Vd||i2(j^) 
By the strong convergence of Vd„ in L'^{Q) and (j4.14p . we conclude that 

Vd„ Vd„ ^ Vd Vd in v' {Q x (0, T)). 

Similarly, 

i|Vd„|2/3 ^ ^iVdpIg in v\Q X (0,T)), 
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aoad 

u„Vd„ uVd in v' {U x (0,r)), 

where we used 

u„ ^ u weakly in L^{[0,T]; HQ{n)). 



Therefore, (|3.ip and (|4.8p hold at least in the sense of distribution. Moreover, by the 
uniform estimates on u, d and (jl.lcp . we know that the map 

t — )• / dn{x,t)ip{x)dx for any ip € ^{0,), 
Jn 

is equi-continuous on [0, T]. By the Ascoli-Arzela Theorem, we know that 

d{x, t)(p{x)dx 



is continuous for any ip € T){^). Thus, d satisfies the initial condition in (j3.ip . 
Now we have the existence of a global solution to (|4.ip as follows: 

Proposition 4.1. Assume that VL is a hounded domain of the class C'^^'^ , v > 0; 

and there exists a constant Cq > 0, such that d • /(d) > for all |d| > Co > 0. Let 
£ > 0,5 > 0, and f3 > max{4, 7} be fixed. Then for any given T > 0, there is a solution 
{p, u, d) to the initial-boundary value problem of (j4.ip in the following sense: 
(1) The density p is a nonnegative function such that 

p(^L\%T\-W^'-{n)), dtpeL-{{0,T)xn), 

for some r > 1, the velocity u G L^([0, T]; i/g (il)), and (14. lap holds almost everywhere 
on (0,T) X Q, and the initial and boundary data on p are satisfied in the sense of traces. 
Moreover, the total mass is conserved, i.e. 

p{x,t)dx = / psfldx, 

for all t G [0,?^]; and the following inequalities hold 

6 [ [ p^+^dxdt < C(e), 
Jo Jn 



I vp'i^ 

Jn 



\\7p\^dxdt < C with C independent of e. 



(2) All quantities appearing in equation (|4.1bp are locally integrable, and the equation 
is satisfied in D (Q x (0,T)). Moreover, 



puGC([0,r];L^,(f])), 



and pu satisfies the initial data. 

(3) All terms in (I4.1cp are locally integrable on 0, x (0, T). The direction d satisfies the 
equation (I3.1ap and the initial data (|3.1bp in the sense of distribution. 
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(4) The energy inequality 

dt ^ Qplup + -l^p^ + ^/ + ^|Vd|2 + AF(d)) dx 

+ 11 f |Vupdx + A / \^d-f{d)\^dx 
Jn Jn 

< 

holds almost everywhere for t E [0,T]. 

To complete our proof of the main theorem, we will take vanishing artificial viscosity 
and vanishing artificial pressure in the following sections. 



5. Vanishing Viscosity Limit 

In this section, we will pass the limit as e — t- in the family of approximate solutions 
(/Oe,Ue,d£) obtained in Section 4. The estimates in Proposition 14. II are independent of n, 
and those estimates are still valid for (jOe, u^, d^). But, we need to remark that will lose 
some regularity when e — )• because the term eAp^ goes away. The space L°°(0, T; L^(Q)) 
is a non-reflexive space, and the artificial pressure is bounded only in space (0, T; (Q)) 
from the estimates of Proposition 14.11 It is crucial to establish the strong compactness of 
the density pe for passing the limits. To this end, we need to obtain better estimates on 
the artificial pressure. 

5.1. Uniform estimates of the density. We first introduce an operator 

S : |/ G LP{n) : l^fdx = 0^^ 

which is a bounded linear operator satisfying 

ll-^[/]llH/,}'''(n) - c{p)\\f\\LP{n) for any 1 < p < oo, (5.1) 

where the function W = B[f] € solves the following equation: 

divl^ = fmn, W\dn = 0. 

Moreover, if the function / can be written in the form / = divg for some g £ L"^ , and 
9 ■ '^lan = 0, then 

\\B[f]\\L-{n) < c{r)\\g\\Lr^n) 
for any 1 < r < oo. We refer the readers to [6l [7] for more background and discussion of 
the operator B. Define the function: 

ip{t,x) =xl:{t)B[p,-p], i;eViO,T), 0<^<1, 

where 

Since p^ is a solution to (|4.1ap . by Proposition 15. II and /3 > 4, we have 

p,-pGC{[0,T],L\n)). 
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Therefore, from dnH), we have G C([0, T], In particular, 93(4, x) G C{[0,T]x 

Q) by the Sobolev embedding theorem. Consequently, ip can be used as a test function 
for (I4.1b|) . After a little bit lengthy but straightforward computation, we obtain: 

Jo Jn 

f-T r r-T r- 

^{pI + 5p^)dxdt + I I iI^PsVl^BIps — ^dxdt 
Jn Jo Jn 

f-T 



+ p I I ipVueV B[pe -p]dxdt 
Jo Jn 

T f 

tpPeUir (g) UeVi3[pe — p]dxdt 



Jn 

T 



e / i'PeUeB[Aps]dxdt (-5 2) 



Jn 

/ 'ippsUsB[dw{psUs)]dxdt 
Jn 

T 



+ e I [ VueVpeB[pe-p\dxdt 
Jo Jn 

~^lo In (^'*' ® ~ ^'^^ + ^(d))^3) ^V5[p, - p]dxdt 



7 

To achieve our lemma below, we need to estimate that the terms Ii — Ij are bounded. We 
can treat the terms related to Pe,Uj similar to [7j. It remains to estimate the term Ij. 
Indeed, 



IM =^11^ (^^^ ® - (^^^ + ^(d))^3) ^VB[pe - p]dxdt 

<C\ \\Vde\\\i(^^-^\\B[pe-p\\\w^'\n)dt + C \\B[pe -p]\\w^'^{n)dt 

<C, 
where we used 

\\B[p£ - p]\\w^'^{n) < CqWps - p\\l^(q), 
Vde £ L\[0,T] X n), 
and /3 > 4. Consequently, we have proved the following result: 

Lemma 5.1. Let {p^,u^,d^) be the solutions of the problem ()4.ip constructed in Propo- 
sition \4^ then 

ll/'e||LT+i((0,r)xO) + II/'£||l/3+i((0,T)xC) < C, 

where C is independent of e. 



(5.3) 
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5.2. The vanishing viscosity Hmit passage. From the previous energy estimates, we 
have 

eAps^O in L2(0,r;VF"^'^(J])) 

and 

eVu,Vpe-^0 in L^{0,T;L^{n)) 

as e -> 0. 

Due to the above estimates so far, we may now assume that 

Pe^pmC{[0,nLl^^,m, (5.4a) 

Us^u weakly in L^{0,T; Wo'^(Jl)), (5.4b) 

p,u, ^ pu in C([0,T],L^,(f])). (5.4c) 

Then we can pass the limits of the terms related to ps , similarly to [7] . It remains 
to show the convergence of d^. Following the same arguments of Section 4, by taking a 
subsequence if necessary, we can assume that: 

d,^d mC{[0,T];Li,,,{n)) (5.5a) 

de^d weakly in L^{0,T; H'^{n)) n L°°{0,T; H^{n)), (5.5b) 

de ^ d strongly in L'^{0,T; H^{n)), (5.5c) 

Vde Vd weakly in L^((0,r) x Q), (5.5d) 

Ade - /(de) ^ Ad - /(d) weakly in L'^{0, T; L^{n)), (5.5e) 

F{de) F{d) strongly in L'^{0,T; H^{n)). (5.5f) 

Consequently, letting e — )• and making use of (j5.4p and (|5.5p , we conclude that the limit 
of (/?e,U£,de) satisfies the following system: 

pt + div(/9u) = 0, (5.6a) 

{pu)t + div(pu (g) u) + VP = ^Au - Adiv(Vd Vd - (^jVdp + F(d))/3), (5.6b) 

dt + u • Vd = Ad - /(d) (5.6c) 



where P = ap] + Sps , here K{x) stands for a weak limit of {i^e}. 

5.3. The strong convergence of the density. We observe that Pe,u^ is a strong solu- 
tion of parabolic equation (|4.1ap . then the renormalized form can be written as 

dtb{ps) + div{b{pe)ue) + ib'{pe)pe - 6(pe))divUj 

(5.7) 

= ediv(xnV6(pe)) - exnb {pe)\Vpe\ 

in D'{{0,T) X R^), with b G C'^[0,oo), b{0) = 0, and b' , b" bounded functions and 
b convex, where xn is the characteristics function of 0,. By the virtue of (jS.Tp and the 
convexity of b, we have 

ip{b' {ps)pe -b{ps)))divu^dxdt < / b{pQ^s)dx + / / dt'4)b{pe)dxdt 
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for any ^ G C°°[0, T], < < 1, V'(O) = 1, ip{T) = 0. Taking b{z) = zlogz gives us 
the following estimate: 

/ / ^/^/Oedivuedxdt < / po,<5 log(po,<5)da; + / / dtil^ log pedxdt, 
Jo Jn Jn Jo Jn 

and letting e — )• yields 



ippdwudxdt < / po^slog po^sdx + / / dfipplog pdxdt, 
Jn Jn ' ' Jo Jn 

that is, 

/ / pdiYudxdt < / po^slog po^sdx — / p\ogp{t)dx. (5-8) 
Jo Jn Jn ' ' Jn 

Meanwhile, {p, u) satisfies 

dtb{p) + div(6(p)u) + - 6(p))divu = 0. (5.9) 

Using (|5.9|) and b{z) = zlogz, we deduce the following inequality: 

/ / pdivudxdt < / po^slog Po,sdx — / plogp{t)dx. (5.10) 
JO Jn Jn ' ' Jn 

From (I5.10|) and ()5.8p . we deduce that 

/ plogp — /5log(/9)(r)(ix < / / pdivu — pdivudxdt (5-11) 
Jn Jo Jn 

for a.e. r S [0, T]. 

To obtain the strong convergence of density ps, the crucial point is to get the weak 
continuity of the viscous pressure, namely: 

Lemma 5.2. Let (pe,Ue) be the sequence of approximate solutions constructed in Propo- 
sition^4jj\ then 

t-T 



lim / / iljr]{ap2 + 5p^ — pdiv\ii;)p^dxdt 
£^0+ Jo Jn 

/ / '4)ri{P — pdlvvi) pdxdt for any ip ^'^{'^^T), ri£'D{Q), 
Jo Jn 



'0 Jn 
where P = ap'i + 5pl^ . 

Proof. We need to introduce a new operator 

Ai = ^-\d^^v), i = 1,2,2,, 

where stands for the inverse of the Laplace operator on M^. To be more specific, Ai 
can be expressed by their Fourier symbol 



M-) = F-\^F{-)), i = 1,2,3, 



with the following properties (see [7]): 

< c(s, 17)||i;||is(K3), 1 < s < oo, 

\\Av\\L9(n)<c{q,s,Vl)\\v\\Ls(pi), q<oo, provided ^ > ^ - ^, 
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and 



Pi^'llL°°{c) < c(s, 0)||w||is(R3) if s > 3. 
Next, we use the quantities 

ip{t,x)=i^{t)r,{x)Mp,], i;eV{0,T), rieV{n), i = 1,2,3, 
as a test function for ()4.1bp to obtain 

/ ipri{{p] + Sp^) - pdivus) pedxdt 
Jo Jn 

= p / ipVusVriA[ps]dxdt - / ip{p] + 6p^)VriA[ps]dxdt 
Jo Jn Jo Jn 



Jn Jo Jn 

T r i-T 



V'/OeUj (g) u^V riA[p^]dxdt — I i ■iptrip^u^A[p^]dxdt 

Jo Jn 

e / ijr]psUsA[diY{xn^ Pe)]dxdt 
Jo Jn 

+ e / / iprjV pi;Vui;A[pi;]dxdt + p / ipUir'VriPsdxdt 

Jo Jn Jo Jn 

-p / iljUsVr]VA[pe]dxdt + / / ipUeipsRlpi.Us] - PeUeR[p£])dxdt 
Jo Jn Jo Jn 

A 



Jn 

T 



^ ^ {vde Vd, - (^|Vd,|2 + F{de))I^ i>Vr,A[p,]dxdt 

£ {Vde Vd, - (^|Vd,|2 + F{de))I^ ^l]VA[pe]dxdt 



A 



where xn is the characteristics function oi Q., A[x\ = V /S. ^[x] 
Meanwhile, we can use 



(5.12) 
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as a test function for ()5.6bp to obtain 
Jo Jn 



^prj{P — iidiYu) pdxdt 
■T r r-T 



[ I i}VvLVriA[p\dxdt - [ [ ipPVrjAlpldxdt 
Jo Jn Jo Jn 

— / ippu^ uVriA[p]dxdt — / / ipt'HP^MPld'Xdt 
Jo Jn Jo Jn 

+ p / i^uVrjpdxdt - p / / ijuVrjVAlpjdxdt (5.13) 
Jo Jn Jo Jn 

r-T 



+ / 'ipu{pR[pu] — puR[p])dxdt 
Jo Jn 

j (^Vd0 Vd- (^|Vd|2 + F(d))l3^ VVr?A[p]dxdi 
^ j [vdQVd-{]^\Vd\'^ + F{d))I^i)r]VA[p]dxdt. 



For the related terms of /Oe, u^, following the same line in [7j we can show that these terms 
in (|5.12p converge to their counterparts in (|5.13p . It remains to handle the terms related 
to dg in (j5.12p . By virtue of the classical Mikhlin multiplier theorem (see [7j ) , we have 



VA[p,] ^ VA[p\ in C{%T]-Li^^^m as e ^ 0, (5.14) 



and 



Since 



A[pe] A[p\ in C((0, T) x Vt)) as e ^ 0. (5.15) 

/ I Vds VdeA[pe] - Vd Vd^[p] | dx 
Jn 

< [ \Vde\'^\A[pe]- A[p]\dx+ [ \Vde\\Vde-Vd\\A[p]\dx (5.16) 
Jn Jn 

+ / \Vd\\Vds-Vd\\A[p]\dx, 
Jn 

using Holder's inequality to ()5.16p . by ()5.14p . ()5.15p . and ()5.5cp we have 

{VdsQVds)'ipV7]A[pe]dxdt^ [ [ {VdQVd)^Vr]A[p]dxdt as e 0. 
I Jo Jn 



T 

Jn 
Similarly, 

/7(- 

Jo Jn 2 



l-\Vde\'^l3)'iljVr]A[pe]dxdt^ [ [ {^\Vd\'^l3)'iljVr]A[p]dxdt as e 0. 
2 Jo Jn ^ 
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Using the strong convergence of F(de), we conclude that, 
A 



^ ^ (vd, Vd, - (^|Vd,|2 + F(d,))/3)) iPVr]A[pe]dxdt 

j (^dQVd- {^\Vd\^ + F{d))h)^^Vr]A[p\dxdt as e ^ 0. 



A 



And similarly, 
A 



^ A 

So we deduce that 



^ ^ (Vde Vd, - (^|Vd,|2 + F(d,))/3)) lPl]VA[pe]dxdt 

j (^dQVd- {]^\Vd\^ + F(d))h)^ijiiSJA[p]dxdt as e ^ 0. 



lim / / ipr]{p2 + 5 — pdivu^) pedxdt 
-^0+ Jo Jn 

= / iprjI^P - pdiYu) pdxdt for any -0 S P(0, T), r/ € 



where P = p^ + 5p^. The proof of Lemma 15.21 is complete. □ 



From Lemma |5.2| we have 



I [ pdivu - pdwudxdt < - [ [ {Pp - ap^+^ + Spl^+^)dxdt. (5.17) 
Jo Jn Jo Jn 

By dsn]) and (f^TTD . we deduce that 

[ plog{p) - p\og{p){T)dx <- [ [ {Pp - ap^+^ + 6pf^+^)dxdt, 
Jn Jo Jn 

and 

Pp - p7+l + < 

due to the convexity of p'^ + Sp^ . So 



plog(p) - plog{p){t)dx < 0. 



On the other hand, 



p\og{p) - plog(/9) > 0. 



Consequently plog(p) = plog{p) that means 

Pe^pin L\{0,T) xn). 

Thus, we can pass to the limit as e — t- to obtain the following result: 

Proposition 5.1. Assume Q, C M.^ is a bounded domain of class C'^~^^,'d > 0. // there 
exists a constant Cq > 0, such that d ■ /(d) > for aZZ |d| > Co > 0. let 6 > 0, and 
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be fixed. Then, for any given T > 0, there exists a finite energy weak solution (p, u, d) of 
the problem: 

pt + div(pu) = 0, (5.18a) 
{pu)t + div(pu (g)u)+ V{p^ + 5p^) 

= iiAu - Adiv (vd Vd - (^|Vdp + F(d))/3^ , (5.18b) 

dt + u • Vd = Ad - /(d) (5.18c) 

with the boundary condition u.\Qfi = 0, d|afi = do and initial condition ()1.2p . Moreover, 
p G L^~^^{{0,T) X O) and the equation (j5.18ap holds in the sense of renormalized solutions 
on D'[{0, T) X M'^) provided p, u were prolonged to be zero on Furthermore, (p, u, d) 

satisfies the following uniform estimates: 

sup \\p{t)\\l-,,^.<CEs[po,mo,do], (5.19) 
te[o,r] ^ ' 

6 sup \\p{t)f, <CEs[po,uio,do], (5.20) 

sup ||v^u(t)||2, <CE5[po,mo,do], (5.21) 

tG[0,T] 

l|u(*)llL2([o,T];//i(n)) < CEs[po,mQ,do], (5.22) 
sup \\Vd\\l2(^) < CEs[po,mo,do], (5.23) 

iG[0,T] 

l|d||L2([o,T];/f2(f7)) < C£;5[po,mo,do], (5.24) 
where C is independent of 5 > and 

Es[po,mo,do]= [ (]-\l^ + J—pl^ + ^pl + hvdo\' + XF{do)) dx. 



Remark 5.1. Recalling the modified initial data (|4.3p - (j4.6p . we conclude that the modified 
energy Es[po,mQ,dQ] is bounded, and consequently the estimates in Proposition 15.11 hold 
independently of 5. 



6. Passing to the Limit in the Artificial Pressure Term 

The objective of this section is to recover the original system by vanishing the artificial 
pressure term. Again in this part the crucial issue is to recover the strong convergence for 
PS in space. 
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6.1. Better estimate of density. Let us begin with a renormalized continuity equation 

b{ps)t + div(6(p5)u5) + {b'{ps)p5 - 6(p5))divu5 = in P'((0, T) x W^) 

for any uniformly bounded function b G C^[0, oo). We can regularize the above equation 
as 

dtSm[b{p)] + div(5„[6(p)]u) + S^[{b'{p)p - 5(p))divu] = g„ on (0,T) x , (6.1) 

where Sm{v) denotes a spatial convolution with a family of regularizing kernels, and 

^ in l2(0, T; as m ^ oo, 

provided b is uniformly bounded (see details in {?]). 

We use the operator B to construct multipliers of the form 

^{t,x)=i;{t)B[Sm[b{p&)]- j^S,n[b{p&)]dxl i;€ViO,T), 0<^P<1, 

where the operator B was defined in Section 5. Taking b{ps) = pg, using (|6.ip and (|5.19|) . 
with fj small enough, we see that 



is in the space C([0, T]; L^(f2)) for any finite p > 1. By (15. ID and the embedding theorem, 
we have (p{t,x) G C([0, T] x Q). Consequently, ^p{t,x) can be used as a test function for 
(|5.18b|) . then one arrives at the following formula: 

[ ^{p] + 5p^s)Sm[Ps]dxdt 
Jn 



T 



1 



iptpsVisBlSmip's] - T7^ Sm[ps]dx]dxdt 
Jn l"l Jn 



+ [ [ ^(/"^Vu^ - (g) U5)VS[S'm[pJ] - / Sm[ps]dx]dxdt 
Jo Jn l"l Jn 

+ / tppsusBlSmips - aps)divu- / S^Kps - (^Ps)'i^'^^5]dx]dxdt 
Jo Jn l"l Jn 

'ippsusB[diySm[{psUs)]]dxdt + [ [ iljpsusB[qm - [ qmdx]dxdt 

Jo Jn l"l Jn 



10 Jn 
+ A 



y^-fi+ / / V'P5U5S[gm - / qmdx]dxdt. 
,_ 1 JO l"l Jn 
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Noting that — ^ in L^(0, T; L^(M^)) as m — )• oo, we can pass to the hmit for m — )• oo 
in the above equahty to get the following: 

6 



Jn 



Now, we can estimate the integrals Ii — Iq as follows. 
(1) We see that 

h 



is bounded uniformly in 5 provided o" < 7 by (j5.19p and (|5.20p . 

(2) As for the second term, by (|5.19p . (j5.2ip . (j5.22p and together with the embedding 

W^^^iVL) ^ L°°{VL) for p > 3, we have 

[■T 

\h\ 



[ iptPs^sBlSmiPs) - iTTi Sm{ps)dx]dxdt 
Jn l"l Jn 

<c [ \ipt\dt < C 
Jo 



provided a < ^. 

(3) Similarly, for the third term, we have 

/ V(/iV'Vu5 - psus Us)VB[S.m{ps) 

) Jn 

< C 



1 

W\ 



Sm{ps)dx]dxdt 



if we choose o" < ^; 

(4) For I4, by Holder inequality, we have 
rT 



where 



/ / ippsusB[Sm{p5 - cr ps )dwu - [ Smips - (TPs)divusdx]dxdt 
Jo Jn l"l Jn 

< / \\P6\\Ly{n)\\^5\\L<^(n)\\Ps^^''^^s\\Li{n)dt, 
Jo 

67 67 



57 — 6 ' 



77-6' 



If we choose cr < ^ — 1, and use (15.191) . (15.201) and (I5.22p . we conclude that 1^ is uniformly 
bounded. 

(5) Using the embedding inequality, we have 



Jn 

T 



i^psusB [div Sm (Ps^s)] dxdt 

< / l|P<5||L^I|u5||L6||j05U5||LP(it 

Jo 

< / IIP^IU^IIu^lliell/oJIlLr-fit, 
Jo 
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where r = . If we choose cr < ^ — 1, and use ()5.19p . ()5.20p and ()5.22p . then Is is 
bounded. 

(6) Finally, we estimate term Ig, let a < ^, then 

= ^V'(vd5 0Vd5-(i|Vd5|2 + F(d5))/3) VS[5„(p^)-^^5™(/5^)dx]dxdt 

<C r \\Vds\\l,^^^\\VB[SM)-^^ I SM)dx\\\L2i^^)dt 

i/ I I "/ 

+ C r \\VB[SM)-^^ j 5n^(p?)d^]||L2(^)d^ 
i/ I I "/ 

<c, 

where we used the smoothness of F, (15. ip . (I5.19p . (I5.20|) and 

Vd^ G L^((0,r) X !^). 
All those above estimates together yield the following lemma: 
Lemma 6.1. Let 7 > |. There exists cr > depending only on 7, such that 

pI'"' + dpi""" is hounded in L\{0,T) x n). 

6.2. The limit passage. By virtue of the estimates in Proposition 15.11 and Remark 15. H 
we can assume that, up to a subsequence if necessary, 

ps^pmCi[0,TiLl^^,m, (6.2) 

us^u weakly in L'^{[0,T]; H^{n)), (6.3) 

d5 ^ d weakly in L^{[0,T]; H^{n)) n L^{[0,T]; H^{n)). (6.4) 

d^ ^ d strongly in L'^{0,T; H^{n)), (6.5) 

Vd5 Vd weakly in L^((0, T) xQ), (6.6) 

Ad5 - f{ds) ^ Ad - /(d) weakly in L'^{0,T; L'^{n)), (6.7) 

F{ds) F{d) strongly in L^{0,T; {Q)). (6.8) 
Letting 5 — )• 0, we have, 

pI^J^ weakly in L^((0,T) x (0)), (6.9) 

subject to a subsequence. 

From (j6.5p and ()6.8p . we have, as 5 — )• 0, 

Vd^ Vd5 - (^|Vd5|2 + F(d5))J3 

1 (6.10) 
^ Vd0 Vd - (-|Vd|2 + F(d))/3 in p'(0 x (0,r)), 

and 

U5-Vd5^u-Vd in P'(17 X (0,r)), (6.11) 

as 5 0. 
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On the other hand, by virtue of (I5.18bl) . (I5l^ - (l5:22| ). we obtain 



psus ^ pu in Ci[0,T];L-2aki^))- (6-12) 
Similarly, we have, as 5 — )• 0, 

d^^din C([0,r];LLafc(f^))- 

By Lemma l6.lt we get 

(5pf ^ in L^((0,r) xQ) asd^O. 

Thus, the hmit of (p, pu, d) satisfies the initial and boundary conditions of (jl.2p and (|1.3p . 

Since 7 > |, (16. Sp and (I6.12P combined with the compactness of H^{Q,) ^ imply, 
as 5 — )• 0, 

(g) U5 pu (g) u in P'((0,T) X ri). 

Consequently, letting 5 — )• in (jS.lSp and making use of (|6.2p - (j6.12p . the limit of {ps, us, ds) 
satisfies the following system: 

Pt + div(pu) = 0, (6.13a) 
{pu)t + div(pu u) + = pAu - Adiv ^Vd Vd - (^jVdj^ + F(d))/3^ (6.13b) 
dt + u-Vd = Ad-/(d) (6.13c) 
in V'{n X (0,T)). 

6.3. The strong convergence of density. In order to complete the proof of Theorem 
12.11 we still need to show the strong convergence of ps in L^{Q), or, equivalently p'y = p"^ . 

Since ps,us is a renormalized solution of the equation (|6.13ap in V ((0, T) x M^), we 
have 

Tk{ps)t + diY{Tk{psus)) + {Tk{ps)p6 - Tk{ps))diY{us) = in v' {{0,T) x R^), 
where Tfc(z) = /cT(f ) for z G M, A; = 1, 2, 3 • • • and T G C°°(M) is chosen so that 

T{z) = z for z < 1, T{z) = 2 for z > 3, T convex. 
Passing to the limit for 5 — )• we deduce that 



dtn{p) + div((rfc(p))u) + (r;(p)p - rfc(p))divu = o in p'((o,r) x m^)), 

where 



Tl{p5)ps - Tk{ps)dvvus {Tk{p)p - Tk{p))d\vu weakly in -L^((0, T) x 9), 

and 

Tkips) ^ n{p) in C{%T]-Ll^^^{n)) for all 1 < p < oo. 
Using the function 

<^{t,x) = ij{t)r]{x)Ai[Tk{ps)l ijeV[0,T], r/ G P(J7), 

as a test function for (|5.18bp . by a similar calculation to the previous sections, we can 
deduce the following result: 
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Lemma 6.2. Let {ps,us) be the sequence of approximate solutions constructed in Propo- 
sition [5A[ then 

lim / / il)ri{pj — fidiYUs)Tk{ps)dxdt = / / il)ri{p'^ — pdwu)Tk{p)dxdt 
s^oJq Jq Jq 

for any S P(0,r),?? G 

In order to get the strong convergence of pg, we need to define the oscihation defect 
measure as follows: 

OSC-,+i[ps ^ p]((0,r) xn) = sup llm sup r [ mips) - T,{pW+^dxdt. 

k>lS^O Jq Jq 

Here we state a lemma about the oscillation defect measure: 
Lemma 6.3. There exists a constant C independent of k such that 

OSC^+i[ps ^ p]{{0,T) xn)<c 

for any k > 1. 

Proof. Following the line of argument presented in [7], and by Lemma 16.21 we obtain 



OSC^+i[ps ^ p]i{0,T) xn)<lim [ [ divusTkips) - divuTfc 



{p)dxdt. 



On the other hand, 



lim / / divu5Tfc(/>5) — divuTfc(p)(ia;dt 

= lim / / [Tkips] - Tk{p) + Tk{p) - Tk{p))divusdxdt 

< 2sup||Vu5||i2((o,T)xQ) limsup||rfc(/)5) - Tfe(p)||i2((o,r)xn)- 

So we can conclude the Lemma. □ 

We are now ready to show the strong convergence of the density. To this end, we 
introduce a sequence of functions G C^(M) : 



Lk{z) = I 

Noting that can be written as 



zlnz, < z < k 

zln{k) + z ^^^ds, z > k. 



Lk{z) = I3kz + bkz, 
where 6^ satisfy (j2.7p . we deduce that 

dtLkips) + div{Lk{ps)us) + Tk{ps)divus = 0, (6.14) 

and 

dtLkip) + div(Lfc(p)u) + rfc(p)divu = (6.15) 
in D'((0,T) X n). Letting 5 — )• 0, we can assume that 

Lkips) ^ L^ip) in Ci[0,T];Ll^^,in)). 
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Taking the difference of (j6.14p and (j6.15p . and integrating with respect to time t, we obtain 
(Lkips) - Lk{p))(l)dx 

n 

= j j (^{Lk{ps)us - Lk{p)u) -^(1)+ {Tk{p)d\vu-Tk{ps)dwus)4>^dxdt, 
for any (j) G Following the line of argument in [7j, we get 



Li 



Lk{p) - Lk{p)j {t)dx 

T r r-T j- (6-16) 

Tk{p)djvvLdxdt — lim / / Tk{ps)d\v\isdxdt. 
•5^0+ Jo Jo. 



Jo. 



We observe that the term Lk{p) — Lk{p) is bounded by its definition. Using Lemma [67 
and the monotonicity of the pressure, we can estimate the right-hand side of ()6.16p : 



Jn 

f-T 



Tk{p)divudxdt — lim / / Ti.{ps)divusdxdt 
s^o+ Jo 



< f f {Tk{p)-Tk{p))divudxdt. 
Jo Jn 



(6.17) 



/o Jn 

By virtue of Lemma 16.31 the right-hand side of (I6.17|) tends to zero as /c — )• oo. So we 
conclude that 



plog(p)(t) = plog(p)(t) 

as A: — )• oo. Thus we obtain the strong convergence of ps in L^((0, T) x Q,). 
Therefore we complete the proof of Theorem 12.11 



7. Large-Time Behavior of Weak Solutions 

The aim of this section is to study the large-time behavior of the finite energy weak 
solutions obtained in Theorem 12.11 

First of all, from Theorem 12. H we have 



esssupi>o^(0 + / / (/u|Vu|2 + A|Ad-/(d)p)dx(it <^(0), (7.1) 
Jo Jn 



where 



m = ^ Qp|u|2 + :f^ + A^lVdp + AF(d)) dx. 

Following the argument in [8], we take a sequence 

Pm{t,x) := p{t + m,x); 

Urn{t,x) := u{t + m,x); 

dm{t, x) := d{t + m, x), 

for all integer m, and t G (0, 1), x G fi. 
From (|7.ip . we have 

Pm G L°"([0,1];L^(J7)), d„, G L°°([0, 1]; (O)), 
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^Vim e L°°([0,l];^'(f^)), PA G L°°([0,l];L^(r?)), 
which are independent of m. Moreover, we have 

lim / (||Vu^||i2(f^) + ||Ad^-/(d^)||i2(f^))(it = 0. (7.2) 

m— i-oo Jq 

So we can assume that, up to a subsequence if necessary, as m — )• oo, 

Pm{t-,x) — >• ps weakly in -^'''((0, 1) x fi); 

vLm{t,x) Ms weakly in L^([0, 1]; i^o (^)); 
Ad„ - /(d„) ^ Ad, - f{ds) weakly in ^^([0, l];L^{n)). 

Furthermore, 

/ Psdx < lim inf / Pmdx < C{Eq). 
Jn Jq 

Using ()7.2p and the Poincare inequality, we have 

lim / llumllla/rndt = 0. 

By the embedding of L^, we have 

= almost everythere in (0, 1) x (7. 

From (|7.2|) again, we have 

Ads — /(ds) = almost everythere in (0, 1) x Q, (7.3) 
under the boundary condition 

d.lsf^ = do. (7.4) 

By the elliptic theory, there exist a unique solution ds £ C^{n) D C(S7) to ([LSI) and (Oil . 
On one hand, from (I7.ip and (I7.2j) . we have 

lim / ( ||pm|Um|^|| _37 + ||Pm|Um||P 67 ] dt = 0. (7.5) 



Since u are the solutions to (jl.lap in the sense of renormalized solutions, one has, in 
particular, 

pt + div(pu) = in V'{{0, T) x 0). (7.6) 

Taking a test function (p{t,x) = ip{t)<j){x) in (|7.6p . where ip{t) G P(0, 1), </> G ^{Q,), we 
have 

/ ( / Pm4>dx I 'lp'{t)dt + / / PmUm'V(l)ljjdxdt = 0. 

Jo \jn / JO 

Letting m — )• cxd and using (17. 5p , we have 

^ (^j^Ps(t>dx^4^'{t)dt = 0, 

which means that ps is independent of time t. 
Similar to Lemma l6.ll we have 

P 



is bounded in L^((0, 1) x $7) independently of m > 0, 
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for some a > 0. So we conclude that 

pln^J^ weakly in ^^((0, 1) x (7.7) 
Therefore, passing to the limit in (jl.lbp and using ()7.2p . (j7.5p . we obtain 

V^=-Adiv(^Vd3 0Vd, -(^|Vd,|2 + F(d,))/3^ \uV'(Sl). (7.8) 

where d^ is the solution to (|7.3p with its boundary condition (j7.4p . 

On the other hand, we can use L^— version of celebrated div-curl lemma to show that 
the convergence in (|7.7p is strong. We refer the readers to [8] and [13] for details. Due to 
the strong convergence in (j7.7p . we have 

Pm Ps strongly in L'^((0, 1) x Q.). 

This, combined with (f7?f|) and fTS]) . gives 

Vp] = -Adiv (vds Vd, - (^|Vd,|2 + F{ds))h 

in the sense of distributions. Denoting 

H = p2- AF(d,), 
and using (17. Sp to rewrite the above equation as follows: 

VH = -AVd,/(d,). 

Notice that 

/(d,) = Vd.m), 

we have 

V(i7 + AF(d,)) = 0, 

that is, 

Vp] = 0. 

Finally, by the energy inequality, the energy converges to a finite constant as t — t- oo: 



:= lim E{t), 



by (|7.5p . we have 

/•m+l 



Thus, 

f>m+l 



lim / / p\u\'^dxdt = 0. 



Eoo= lim / / {-p\u\^ + -p^ + A-|Vd|2 + AF(d) 

= X (^''^ + ^^iVd.P + AF(d,)^ dx. 

lim / / F{d)dxdt = / F{ds)dx 



dxdt 



We observe that 

r>m+l 



and 



m+l 



lim / / -\\7d\^dxdt= / -|Vdspdx 
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because of (|7.3p . Moreover, using (jl.lap one can easily see that 

p(t,x) — )■ ps weakly in L'^(il) as t ^ oo. 

Thus, we have 



which means 



< lim inf [ f />^ + A-|Vdp + AF(d) ) da; 

7n V7 - 1 2 J 

< lim sup / (^—p2 + A^|Vd,|2 + AF(d,)) dx 

Ih^T / f ip|u|2 + -^pj + Ai|Vd,|2 + AF(d,)) dx 
^ooj^ \2 7-12 J 



< 

t 



lim = 



hm / -^—p'^dx = / ^—p^dx, 

2.8p follows since the space L"^ is uniformly convex. 
The proof of Theorem 12.21 is complete. 
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